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The recent discovery of topological Kondo insulators has triggered renewed interest in the well-known Kondo 
insulator samarium hexaboride, which is hypothesized to belong to this family. In this Letter, we study the spin 
texture of the topologically protected surface states in such a topological Kondo insulator. In particular, we 
derive close relationships between (i) the form of the hybridization matrix at certain high-symmetry points, (ii) 
the mirror Chem numbers of the system, and (iii) the observable spin texture of the topological surface states. In 
this way, a robust classification of topological Kondo insulators and their surface-state spin texture is achieved. 
We underpin our findings with numerical calculations of several simplified and realistic models for systems like 
samarium hexaboride. 

PACS numbers: 71.27.+a, 75.20.Hr, 73.20.At, 03.65.Vf 


Introduction — Since the theoretical characterization of 
topological Kondo insulators (TKls) [1, 2], this class of 
materials has attracted much attention in the community. 
One material in particular, samarium hexaboride (SmBg), 
has been studied extensively both theoretically and exper¬ 
imentally. Several theoretical studies predicted SmBg [3- 
8] and related compounds [9, 10] to be TKIs with pro¬ 
tected gapless surface modes. Different experiments showed 
that, at sufficiently small temperatures, transport is in¬ 
deed dominated by the surface contributions [11-13]. At 
the same time, angle-resolved-photoemission-spectroscopy 
(ARPES) [14—19], quantum-oscillation [20], and scanning- 
tunneling-microscopy measurements [21] confirmed the exis¬ 
tence of gapless surface states. Nevertheless, due to the small 
bulk gap of 15-20 meV [22-24] and strong electronic corre¬ 
lations, a detailed characterization of the nature of the surface 
states is difficult and may require additional concepts such as 
atomic reconstruction [25], Kondo breakdown [26], or exci- 
tonic scattering [27]. Some groups also challenged the sce¬ 
nario of a TKI [24, 28]. To date, the most conclusive evidence 
for the topological nature of the surface states is provided by 
spin-resolved ARPES measurements of the (001) surface [23] 
showing that the surface states around the X point of the sur¬ 
face Brillouin zone (sBZ) are spin-polarized. 

SmBg is predicted to have a band inversion at the X high- 
symmetry points (HSPs) [3-6]. The X-inverted phase has a 
nontrivial strong Z 2 index vq = 1, weak topological indices 
V = (1,1,1) and protected surface Dirac cones as shown 
in Fig. 1 (a) for the (001) surface. The experimental work 
in Ref. 23 is consistent with these predictions and further¬ 
more suggests that the spin texture of the surface states is as 
sketched in Fig. 1 (c). Interestingly, however, several theoreti¬ 
cal studies reached conflicting conclusions about the nature of 
the spin texture [30-32], which is not uniquely determined by 
the Z 2 invariants. In fact, for linear Dirac cones, two situations 
are compatible with the cubic symmetry, see Figs. 1 (c) and 
1 (d). They are distinguished by opposite winding numbers 

Wx = ±1 of the planar unit spin (n_t, ny) - {SSy)l + Sy 
around the X point, where the winding number around the 
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FIG. 1. (a) Sketch of the Dirac cones in the (001) surface Brillouin 
zone, (b) Positive directions of mirror invariant lines in the sBZ for 
the (001) surface with outward pointing normal vector = e,, see 
Ref. 29 for further details, (c-d) Sketch of the spin (or pseudospin, 
see page 2) textures in the (001) sBZ. While at E the winding number 
is always Wf = 1, at the X points it can be Wx = +1 (c) or vrx = -1 
(d), depending on the configuration of the MCNs. 


HSP K [33] is defined as 

= 2 ^ V [im log(n^ H- iny)] ■ ds, 


Wk 


( 1 ) 


with jK a contour encircling K in an anticlockwise fashion. 
This discrepancy between different theoretical models and ap¬ 
proaches raises the important question of what determines the 
spin texture in cubic TKIs. 

In this Letter, we provide two answers to this question: 
First, we show that there is a close connection between the 
spin texture and the mirror Chem numbers (MCNs) [34]. In 
particular, knowledge of the MCNs allows us to distinguish 
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FIG. 2. Chiral surface states with mirror eigenvalues ±i along the FX 
(a) and XM line (b) in positive direction (see Fig. 1) for Co = -2 and 
C,= 1. 

between the two situations shown in Figs. 1 (c) and 1 (d). Sec¬ 
ond, we provide analytical expressions relating the surface- 
state spin texture to the hybridization parameters of specific 
models. These relations demonstrate that the number and type 
of included orbitals in the effective model does not uniquely 
define the winding number; instead, the relative strength of 
different-range hybridization parameters is equally important. 
In addition, we show how the system can be tuned across topo¬ 
logical phase transitions, during which the surface-state spin 
texture changes while all the Z 2 invariants remain unaffected. 

In the remainder of the Letter, we will provide the details to 
the above statements. We will also apply the general argumen¬ 
tation to a multiorbital model with itinerant and localized 
Fg electrons, as in Ref. 32. Other models are discussed in the 
Supplemental Material [29]. 

Mirror Chern numbers define pseudospin texture — To 
start, we review certain facts about the MCNs in SmBe. The 
MCNs are topological invariants, which are protected by mir¬ 
ror symmetries [8, 34, 35]. In a cubic system, there are three 
distinct MCNs: Co = ^ ^ 

with O',/? 6 {x,y,z} and j3 a, where refers to the 
Chern number of the Bloch states on the mirror-invariant 
plane S with eigenvalue -Hi under the mirror operation, see 
also Ref. 35. As was shown in Ref 8, the cubic symmetry im¬ 
plies that the MCNs in the X-inverted phase are Co = 2 mod 4, 
C;r = 1 mod 4 and Cd = 1 mod 2. These values imply two 
additional Dirac nodes along the FX line on the (110) sur¬ 
face [8, 35]. In the following, we show that the MCNs also 
determine the spin texture on the (001) surface. (A related ar¬ 
gument for Hg-based topological insulators was presented in 
Ref 36.) 

The projections of the mirror planes onto the (001) sur¬ 
face correspond to the high-symmetry lines (HSLs) shown 
in Fig. 1 (b). Along these mirror invariant lines (MILs), 
we can classify the surface states according to their mirror- 
eigenvalues +i. The bulk-edge correspondence for each 
mirror-invariant plane then states that the MCN C is equal to 
the number of right-moving (C > 0) or left-moving (C < 0) 
surface modes with mirror-eigenvalue -Hi, see Fig. 2. There 
exists a certain freedom to choose signs in the calculation of 
the MCNs. We use a convention [29], which leads to the posi¬ 
tive directions shown in Fig. 1 (b). 

The mirror eigenvalues also define a pseudospin of the sur¬ 


face states p in the following way: On the ky = 0 or ky - n 
MIL, we choose a basis jHi,H 2 ) in which the mirror operator 
takes the form M,. = -tpy, where pa is the cr-th Pauli matrix. 
Furthermore, on the - 0 and k^ - ti MILs we can choose 
the mirror operator - -ipx- The pseudospin is then given 
by the spinor u — au\ + bu 2 = {a,by. Its relation to the 
physical spin of the electron is detailed on page 4. It follows 
that, along the MILs, the pseudospin lies in the surface plane 
and is always perpendicular to the MIL. In order to make the 
connection to the pseudospin texture, it is useful to consider 
the effective Hamiltonian close to the Dirac node at the HSP 
X = F or X = X: 

HKiq) = v'KPyqx - v^icPxqy = i{v]^Myqx - v^Mxqy). ( 2 ) 

Here, we measure the momentum relative to the respective 
HSP, q - k - K. At the f-point, the cubic symmetry implies 
that Vp = and the resulting pseudospin texture necessarily 
has a winding number Wf = 1. But at the X points, v;! + v- 
in general, and the winding number of the pseudospin texture 
is Wx = sgn(v^vA). Because the MCNs fix the direction of 
the pseudospin at the points where the Fermi lines cross the 
MILS, the MCNs also fix the relative sign between and vL 
and hence the winding number Wx- It is then easy to see that 
the set (Co, C„) — (2,1) implies the pseudospin texture shown 
in Fig. 1 (c), while (Co, C„) = (-2,1) implies the pseudospin 
texture shown in Fig. 1 (d). For linear Dirac cones at f and X, 
there are no other possibilities, i.e. higher MCNs imply addi¬ 
tional Dirac nodes along HSLs [29]. (Note that in Fig. 1 we 
assume a chemical potential above the Dirac nodes.) 

Hybridization matrix defines mirror Chern numbers — We 
now analyze the connection between microscopic parameters 
of the electronic Hamiltonian and the set of MCNs. From ab- 
initio calculations [5, 30, 32, 37] it is known that the states 
near the Fermi energy in SmBe are predominantly formed by 
the Sm 5d electrons of Eg symmetry and the Sm 4/ electrons 
in the J - 512 multiplet. The latter splits further into a Fg 

quartet, |F<|^) = |+f) + ^|+|) and |F® ) = |+i), and 

a Fy doublet, |r 7 _±) = where the index 

+ is the orbital pseudospin [38]. Our strategy is to start in 
the trivial insulating phase without band inversion and con¬ 
sider the effective model, which describes the gap closing and 
subsequent band inversion at the X points. 

The little co-group at the X point is isomorphic to the tetrag¬ 
onal symmetry group D 4 h. Thus, all the irreducible represen¬ 
tations are at most two-dimensional and the band inversion 
occurs between the energetically highest single Kramers pair 
of / electrons fx,± and the energetically lowest single Kramers 
pair of d electrons t/x.Ti- Near the transition between the trivial 
and the topological phase, the low-energy electronic structure 
can be obtained from an effective 4x4 Bloch Hamiltonian 
around the X points. 
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Equation (3) is given for a spinor if/ = (c/x,T>'^x,i-/x,+ ,/x,-y 
and q is measured from X. The simultaneous presence of 
inversion and time-reversal symmetry allows us to choose 
the hybridization matrix in the form <1)^ = i0^ ■ cr, with 
- 4>q - -<l> q and O' the Pauli matrices in spin space. In 
the following, we consider X = (0,0, tt), and expand to lowest 
order in = e/1 and 

d), = i + o-yqy) + (f>20-zqz\ ■ (4) 

As we show below, the relative sign between the two inde¬ 
pendent parameters and (p 2 of the linearized hybridization 
matrix (4) determines the set of MCNs and hence the surface- 
state spin texture in the X-inverted phase. 

First, we address the MCN Co and therefore consider the 
mirror plane kx -0. The mirror operator in the basis of Eq. (3) 
is Mx = -iTj ® (Tx- Thus, in the subspace Mx - -i-i, Eq. (3) 
reduces to 

= M - (piqyflx + <p2qzMy “ AyUj , (5) 

where jda are the Pauli matrices acting on the basis vectors 
(1, -1,0,0)/ V2 and (0,0,1,1)/ V2 and we have defined s = 
j(sd+Sf) and A s i(e/-e^). The total Berry flux contribution 
of the lower band of a Dirac model h(k) = e(k) d{k) ■ fi with 
d - d/\d\ is 

= , 6 ) 

which for our case with d(k) - {-(I)\ki,(p 2 k 2 , -A) leads to 

C°‘™ = i sgn(A^^i02) ■ (7) 


Therefore, starting from the trivial phase with A < 0 and 
creating a band inversion at X (A > 0) leads to a MCN of 
Co = 2 sgn( 0 i 02 ), where the factor 2 comes from the fact that 
there are two X points in the kx -0 plane. 

The two other MCNs can be calculated analogously, see 
Ref. 29 for details. We obtain C;r = 1 and Cd = vsgn((^i^ 2 ), 
where v = -1 for a band inversion between (x^ - y^) and a 
linear superposition of Eg*^ and r 7 , and v = 1 for a band in¬ 
version between (3z^ - r^) and Eg^^ Hence, if sgn((;ii02) = 1 
(-1), we recover the set of MCNs which imply the pseudospin 
texture in Fig. 1 (c) [Fig. 1 (d)]. In general, we obtain 


v^x = sgn(4>i(f>2). 


( 8 ) 


Model calculations for SmBf, — In the following, we will 
illustrate our theoretical findings by calculations with an ef¬ 
fective lattice model for SmBe. In the interest of simplic¬ 
ity, we will restrict ourselves to the Eg quartet for / electrons 
and study a model similar to that used in Ref. 3 and Ref. 32. 
Analogous calculations can be performed for the full or the r 7 
model [29]. The Bloch Hamiltonian is an 8 x 8 matrix 


77(r8) - 


(hd 

1<1>8 



(9a) 


where the hopping of d and / electrons and the hybridization 
are given by 


hd(k) = 0-0 

hi(k) - 0-0 
08(A:) = -i 






|(C1 + C2)(tf> + 2t'fc3 ^(Ci - C2) tf' - 2tf’c3 

C2) - 2tfc3) - 4 tf C.C2 - 2 t<')c 3 - i(c. . C2) {t^J 

eg - |(ci H- C2) ( 4 '' -H 24^'C3) ^(ci - C2) ( 4 '^ - 24^^C3) 

^(ci - C2) ( 4 '^ - 24^^C3) eg - 44^^0102 - 24^^03 -\{ci+ C2) ( 4 '' 

3/2 Vg^^SlO-l -I- S2O-2) -H 3 y® [(Ci H- C 2 )S 30'3 C 3 (iiO-i H- S 20 ' 2 )] • ■ ■ 

- V3/2 Vg'^sio-i - *20-2) -I- [(ci - 02)530-3 -I- 03(510-1 - 520-2)] ... 

- V 3 / 2 y 445 io-i - 520-2) -H V 3 y® [(oi - 02)530-3 h- 03(510-1 - 520-2)] 
y*'* [2530-3 -H 1/2(510-1 -H 520-2)] -t- y® [(01 -H 02)530-3 + 4(025iO-i + Oi530-3) -H 03(510-1 -H 530-3)] 


(9b) 

(9c) 

(9d) 


with the Pauli matrices o-^ acting in spin space and the spinor 
f = (c/^2_3,2^ll,(i3j2_r2,tj.,/r(i)_j,,/r(2) ^)‘. Here, Ca = cos ka and 
5 q, s sinA:Q,, and we use (yd4)) to denote first and 
second neighbor hopping (hybridization) parameters, respec¬ 
tively. The hopping and hybridization parameters should be 
considered as renormalized due to a strong local Coulomb in¬ 


teraction for the / electrons [39^1]. As long as the elec¬ 
tronic states near the Fermi energy are well described by 
quasiparticles, the adopted single-particle approach to com¬ 
pute the topological invariants is justified, even in the presence 
of strong electron correlations [42^4]. A typical band struc¬ 
ture in the X-inverted phase (without hybridization) is shown 
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FIG. 3. Band structure without hybridization (a) and phase diagram 


-0.03, = 0.02, and eg = -3. The two spin textures (c-d) in phases 

I and II, respectively, are realized for the hybridization parameters 
Vf) = (-0.1,0.!) and (F'b, V®) = (0.3,0.07), respectively. 


(b) for the Fg model defined in Eq. (9) with ij'’ = 1, = -0.2, t^’ = 


,(i) _ 


in Fig. 3 (a). 

For kx - ±ky, the ofFdiagonal elements of both and /tg 
vanish and the d and / electrons are split into (x^ - y^) and 
(3z^ - r^) orbitals, and Fg*^ and Fg^\ respectively. Therefore, 
at the point X = (0,0, n) we obtain 

hdX) = CTO diag [-3 (4'^ - 2tf ), 4') - 2tf] , (10) 

and similarly for the Fg orbitals. Ab-initio calculations [9] 
suggest that > 0 and ^ band 

inversion occurs between the and the Fg ' orbitals. The 

hybridization matrix for these two orbitals can be expanded to 
first order at the X point; 

-i o, = + a-yqy) - 2vf^')-6vf’ a-^q^. (11) 

Therefore, according to Eq. (8), we obtain 

wx = -sgn[yf (y''^-2yf)], (i2) 

leading to the phase diagram shown in Fig. 3 (b). As dis¬ 
cussed above, v = -1 for (x^ - y^) and Fg*^ orbitals, such 
that we expect (Co,C;r, Q) = (2,1,-1) in phase I, lead¬ 
ing to a pseudospin texture with Wx = 1, while we expect 
(Co,C;r,Cd) = (-2,1,1) and Wx = -1 in phase II. At the 
phase transitions Vg =0, the hybridization vanishes along 
the FX line, for V® = | Vg^^ it vanishes at both the XM and 
XR lines. This causes the hybridization gap to close and the 
MCNs (Co, C„, Cd) to change by (±4,0, +2). We numerically 


confirmed the phase diagram in Fig. 3 (b) by directly calculat¬ 
ing the MCNs using a method for a discretized BZ [45]. Fig¬ 
ures 3 (c) and (d) show the physical-spin texture in phases I 
and II, respectively. They were calculated for a slab of 500 
unit cells and fit the expected texture for the pseudospin. 

Relation between physical spin and pseudospin — The 
observed equivalence between physical-spin and pseudospin 
texture in Fig. 3 requires more attention; Because the / elec¬ 
trons experience strong spin-orbit coupling, the orbital pseu¬ 
dospin defined above is not equivalent to the physical spin of 
the electrons and the mirror and spin operators do not com¬ 
mute. The relation between physical and orbital pseudospin 
for the 7 = 5/2 multiplet is given in Ref 29. 

According to the definition of the pseudospin above, a 
surface pseudospin in positive n direction corresponds to an 
eigenvalue -i of M„. In order to find a relation between the 
physical-spin and pseudospin texture of the surface states, we 
therefore consider the effect of the projector s f (1 h- iM„) 
on the physical-spin operator 5„, where projects onto the 
subspace M„ = -i and n is the normal vector of the mirror 
plane. One can show that, for the and J - 5/2 multiplets, 

= 0 for n±n', (13) 

which states that on a MIL, the physical spin is always parallel 
(or antiparallel) to the surface-state pseudospin. Whether the 
two are parallel or antiparallel is determined by the eigenval¬ 
ues of the projected spin operator, 

= p^:snP^:. (u) 

Eor the d orbitals we have S - cr leading to eigenvalues 
- 1-1 of while for the F 7 , the Fg, and the full model, we 
obtain the (approximate) spectra (-0.24), (0.52,0.14), and 
[0.71,0.14,-0.43], respectively, see Ref. 29. As all eigenval¬ 
ues are positive for the Fg model, the physical spin is indeed 
always parallel to the surface-state pseudospin and all find¬ 
ings concerning the pseudospin are directly transferable to the 
physical spin. This is not the case if we also consider the F 7 
orbital, because the projected spin operator of / electrons also 
has negative eigenvalues. In these cases, the relation between 
pseudospin and physical spin of the surface states depends on 
the orbital character of the state. In all cases we have studied, 
the winding number of the physical spin sufficiently close to 
the Dirac node is nevertheless identical to the winding num¬ 
ber of the pseudospin. However, the direction may be reversed 
around some of the Dirac points. Indeed, we find that this may 
occur for the F 7 model, signaling a dominant (in terms of spin) 
F 7 character of the surface states [29]. 

Finally, we mention that for other models with band cross¬ 
ings along some HSLs, there is the possibility of phases with 
higher MCNs and a larger number of protected surface states. 
We discuss an example in Ref. 29. 

Conclusion — We have derived a close relationship be¬ 
tween the hybridization matrix at the X high-symmetry points, 
the mirror Chern numbers, and the spin texture of the topolog¬ 
ically protected surface states in topological Kondo insulators. 
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Although we have motivated our study with SmBg, the line of 
argumentation also applies to other topological insulators. Ex¬ 
plicit calculations for different models for SmBe showed that 
the spin texture of the surface states does not only depend on 
the orbitals that are included in the effective model, but also 
depend on the magnitude of different hybridization parame¬ 
ters. This fact needs to be kept in mind when interpreting 
ab-initio or effective-model-based calculations for this type 
of materials. Finally, our results can be used to infer mirror 
Chern numbers from spin-resolved ARPES measurements and 
predict further observables. 

Acknowledgments — We would like to thank T. Neupert, 
M. Shi, and N. Xu for inspiring discussions. This work is 
financially supported by a grant of the Swiss National Science 
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Note added — During the submission process, a related 
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Supplemental Material for “Surface-state spin textures and mirror Chern numbers in topological 

Kondo insulators” 


MIRROR OPERATORS AND SIGN CHOICE OF THE 
MIRROR CHERN NUMBERS 


In general, the mirror operator for a plane with normal vec¬ 
tor n can be written as 


M„ = IRr, = IRfR^: , 


(SI) 


where I is the inversion operator and denotes a rotation by 

n around n. The spin part of the rotation is given by R'^ - 

-i(n ■ (t), where cr are the Pauli matrices acting in spin space. 

f?™*’ denotes the orbital part of the rotation and depends on 

the symmetry of the considered orbitals. For example, for the 

orbitals > fri > /r(>)> the orbital part is = 

\ if n - Ca {a - x,y,z); but for the mirror planes — ±ki} 

this operator contains additional nontrivial signs. Specifically, 

for n - -^(Cy - Cx) we obtain - viif/i with v = 1 

in the subspaces and v = -1 in the subspace 

[dx2-yl,fym,fY^]■ 

8 

We note that the signs of the MCN are not uniquely deter¬ 
mined in general: One can choose the sign of the mirror oper¬ 
ator M and the orientation of the mirror invariant plane for the 
calculation of the MCN, n^p. We fix the signs of the MCNs by 
the convention 


--i(nmp-tr). (S2) 


The choice of orientation iimp and the sign convention dur¬ 
ing the calculation of the Chern numbers still affects the pos¬ 
itive directions in the sBZ. We use the convention where the 
Berry connection. Berry curvature, and Chern number are de¬ 
fined as - i{t^{k)\dp\il/{k)), Tpv - dpJly - dyJlp, and 
C - ^ f Aki Ak 2 ’F\ 2 , respectively. Note that this convention 
leads to an additional factor of -1 compared to Ref. S1 which 
we use for our numerical calculations of the MCNs. 

For the given conventions for the definition of the mirror 
operator and the Chern number, the positive directions in the 
sBZ can be defined as Mpos = «sf x «mp^ where is the out¬ 
ward pointing normal vector of the surface. In this work, we 
use the conventions n^p = for a mirror-invariant plane 
ka -Qo^ka and n^p = ~ for the plane kx - ky. 

This leads to the positive directions on the top (n^f - +e^) 
(001) surface as shown in Fig. 1 (b). 


MIRROR CHERN NUMBER IN THE X-INVERTED PHASE 

In the Main Text we have calculated the MCN Cq from an 
expansion of the hybridization matrix around the X points. In 
the following we will explicitly show the analogous calcula¬ 
tions for the two other MCNs. 

For the mirror plane k^ - tt, the mirror operator is = 
-ixj ® CTj and the effective Hamiltonian in the - -Hi sub¬ 
space is 

^effl=;r(^) = fil -H (pMyHx + qxRy) “ AjUj , (S3) 

using the basis vectors (0,1,0,0) and (0,0,1,0). It only de¬ 
pends on the parameter (pi and the contribution to the total 
Berry flux amounts to 

^ sgn(A cp\) = i sgn(A). (S4) 

Using the same argument as for Cq and with the fact that there 
is only one X point in the kx - ti plane, the MCN is therefore 
always = 1 in the X-inverted phase. 

Finally, consider the mirror plane kx - ky with Mxy - 
-ivTj - o'x), where v = +1 is the orbital rotation eigen¬ 

value in the considered subspace. There, we obtain the effec¬ 
tive Hamiltonian 

= Sl -H Vpiqxyfix + (piqzRy - (S5) 

where we use qxy ■- ^'fiqx with qx — qy and we chose the 
basis vectors ((1 -Hi)/ 2 , v/ V2,0,0) and (0,0, (1 -Hi)/ 2 , -v/ V2). 
The choice n^p = corresponds to ki = q^ and 

k 2 = qxy- Then, analogous to Eq. (7), we obtain the Berry flux 

= ^^sgn(A0i^^2), (S6) 

corresponding to a MCN Cq - vsgn(0i(^2) in the X-inverted 
phase. 


SIMPLE MODEL WITH NNN HYBRIDIZATION 


In Ref. S2, we have defined a simplified two-orbital model 
to describe SmBg. In order to be able to discuss different spin 
textures of the surface states, we add a NNN hybridization 
term to it. Then, using the spinor ip - the 

Bloch Hamiltonian is given by 




'hf{k) 0(*>(A:)' 

¥^\k) hf(k)^ 


(S7a) 


where hopping and hybridization are defined as 



S2 


h^d^k) = \-2td (ci +C2+ C3) - Atj (ci2 + C 23 + C31) - St'J C1C2C3] 1, (S7b) 

hj^(k) = [e/ - 2 f/(ci +C 2 + C3) - 4f^(ci2 + C23 + C31) - 8f^'ciC2C3] 1 , (S7c) 

= -2 [cr;,S;j. (Vi + y2(c3, + cj) + CTySy (Vi + V 2 (c^ + cS) + o-^s^ (Vi + y2(C;i- + c,,))] , (S7d) 


TABLE SI. The two possible configurations (up to an overall minus 
sign) for MCNs in the simple model with only NN hybridization. 
Multiple band inversions are additive for the MCNs. The first version 
has already been discussed in Ref. S2. 


band inversion 

(Co, Cj., Cd)i 

(Co, C^, Cd )2 

T 

(1,0,1) 

(1,0,1) 

X 

(-2,1,-!) 

(2,1,1) 

M 

(l,-2,-l) 

(1,2,1) 

R 

(0,1,1) 

(0,1,1) 


TABLE SII. Values for Vi and V 2 and conditions for k such that 
OfA) = 0 according to Eq. (S7d) for the simple model (HSPs are not 
shown). A* is a free parameter and permutations of the A„ are always 
allowed. Also shown are the HSLs to which the shown values for A 
correspond. 


Vi 

A 

HSL 

0 

(0, K, A*) 

XM 

-2y2 

(0,0,A*) 

TX 

+ 2 V 2 

(tT, 7T, a*) 

MR 

-Vid+cosk") 

(0,A*,±A*) 

TM 

+ 1 ^ 2(1 - cos A*) 

(n. A*, ±A*) 

XR 

-2y2 cos A* 

A 6 {±A*p 

TR 


with the definitions Ca = cos(ka), Cap = CaCp, and 
Sa = sin(AQ,). The band structure of this model without hy¬ 
bridization is shown in Fig. SI (a). In this simple model, the 
rotation operator is R°* - 1 for all mirror planes. If there 
is only a NN hybridization (Vj = 0), then 0(A) = 0 is only 
possible at one of the HSPs, A e {0, ; 7 r)^. Therefore, only two 
possible sets of MCNs are possible, see Table SI. 

However, in general there are multiple other possibilities 
for 0(A) = 0, if we allow for nonzero NNN hybridization, see 
Table SII. If there exists a band crossing of the bare bands 
on one of the HSLs, phase transitions are possible by chang¬ 
ing the model parameters Vi and y 2 . While the Z 2 topolog¬ 
ical indices are invariant as no band inversion at the HSPs is 
changed, the MCNs can change with this procedure. This leads 
to a richer phase diagram which includes phases with higher 
values for the MCNs and a large number of protected gapless 
surface modes. An example with two additional phases for 
fixed hopping parameters is shown in Fig. SI (b). 


The phase III exists due to the fact, that the gaps at the XM 
and XR lines do not close simultaneously as in the F 7 and Fg 
models. The line separating it from phase II is obtained from 
calculating the crossing point of d- and /-electron bands along 
the XR line and using the formula from Table SII. Similarly, 
the lines separating phases I and III are obtained by calculat¬ 
ing the two crossings of the d and / bands along the FM lines. 
The two other lines, V 2 - -jVi and Vi = 0 correspond to gap 
closings at the FX and XM lines, respectively, similar to the 
results from the Fg model in the Main Text. The MCNs for all 
different phases were calculated numerically using the same 
method as in the Main Text [SI]. 


(b) V 2 



FIG. SI. (a) Bandstructure without hybridization for the model de¬ 
fined in Eq. (S7) with G = F - “0.5, t'J = 0, If = —0.2td, and 
£/ = -2. The hand crossing along the TM line combined with the 
fact that the hybridization can vanish along all HSLs leads to a rich 
phase diagram (h). Phases I and II are the same as for the Tg model 
in the Main Text. For the two additional phases III and IV, the MCNs 
(Co, C„, Cg) are (-2, -3,1) and (-2,1,3), respectively. Note that all 
phases are consistent with the general results of Ref. S3. 


REALISTIC MODELS FOR SMBo 

As discussed in the Main Text, for a more realistic descrip¬ 
tion of SmBg, usually the two Eg orbitals and one or several 
of the spin-orbit coupled / orbitals F 7 and Fg are used [S4- 
S7]. The complete Hamiltonians can be constructed from the 
individual intra- and inter-orbital hopping amplitudes. Here, 
we start from the model defined in Ref. S7 with a selection of 
nonzero parameters as a basis for our numerical calculations. 
In the following definitions, Tq. and ctq, will denote the Pauli 
matrices in orbital and spin space, respectively. Also we will 
again use the definitions Co = cos Aq. and Sa = sin Aq,. 

Then, with hopping and hybridization defined by 











S3 


hd{k) = o-Q 
hi{k) = ctq [e? 
hik) = cTo 


-|(ci + C2) + 2 f®C 3 ) ^(ci - C2) ( 4 ^^ - 24^'C3) 

V^(ci - C2) ( 4 ‘’ - -AtfciC2 - 2 tfc 3 - |(ci + C2) ( 4 ^^ + 24^’c3) 

- 24 ^Vci +C 2 + C3) - 44 ^ 4 ciC 2 + C2C3 + C3C1) - 84^'ciC2C3] , 

1 - i(ci+C 2 ) ( 4 '^ + 24 ^^ C 3 ) 4 ^(ci - C 2 ) ( 4 '' - 24 ^'C 3 ) 

^(Ci - C2) ( 4 '^ - 24^^C3) eg - 44 ^'ciC 2 - 24'^C3 - j(ci + C2) ( 4 '^ + 24 ^^<^3) 


hisik) = ( 4 * 8^0 (-2c 3 + Cl + C2) + 2 tf^ ((ci + C2)C3 + - 2 ciC 2 ) ctq + 2 S 2 (sicr 2 - S 20 'l) 

V 348 Vo(ci - C2) - (2 V 3 cro(ci - C 2 )C 3 - 4 icr 3 SiS 2 + 2 iS 3 (cr 2 Si + CriS2))) , 

<i’7(k) = -i (2530-3 - (sicri + 520-2)) + (2 ^/^CsisiCTi + 520-2) - 2 V 3 (ci + 02)530-3) ... 

y)''(-V 3 ( 5 io-i - 520-2)) + yf 4 c 3 (sio-i - 520-2) -i- 4 (c 25 io-i -C1520-2) - 2 (ci -02)530-3)) , 
3/2 y^Vsio-i + 520-2) + 3 y® [(ci + 02)530-3 + 03(510-1 + 520-2)] ... 

• V3/2 yg^ 45 io-i - 520-2) -H V 3 y® [(ci - 02)530-3 + 03(510-1 - 520-2)] ... 

- V3/2y445io-i - 520 - 2 ) + V3y® [(ci - 02 ) 530-3 + 03 ( 510-1 - 520 - 2 )] 
y^' [2530-3 + 1/2(510-1 + 520-2)] H- y® [(ci + 02)530-3 + 4(C25 iO-i + Ci 530-3) -H 03(510-1 + 530-3)] 


<t>&(k) = -i 


(S 8 a) 

(S 8 b) 

(S 8 c) 

(S 8 d) 

(S 8 e) 

(S 8 f) 


and the spinor if/ = (c/^ 2 _j, 2 _||,c/ 3 j 2 _r 2 the 

Hamiltonian of the full model can be written as 


^full - 


7 dy Ay) _ ~ fy Ay) _ 

’ I7 — ’ ^8 — 


'hd 

(hg^ 

(P 7 I' 



<E)g 

/zg 

hy^'' 


(S9) 


hyi 

hy ^ 



parameters are (for y - 

l,2,3)4t'^- 


(r) 

78 “ 




- vrilJ, ey - and we set e^ = 0 . 


Ft 


Reduced models including only the r 7 or only the r 8 or¬ 
bitals can be obtained by removing one line and column of 
the Bloch matrix in Eq. (S9). 


CALCULATIONS FOR THE Ty MODEL 

Now we want to consider only the r 7 doublet for / electrons 
and the quartet for d electrons. Then, the Bloch Hamilto¬ 
nian is given by 


Kry) = 


hd 

O 7 


07 '' 
hy ^ 


(SIO) 


where the hopping and hybridization parts are defined in 
Eq. (S 8 ). 

As in the Main Text, we consider the situation where the 
(x^ — y^) orbitals are lower in energy at the point X = (0,0, n) 
such that the inversion occurs between those and the r 7 or¬ 
bitals. The hybridization matrix for these two orbitals can be 


expanded to first order at the X point: 

-i % = + (r,qy) (y)'^ + 2 V3y®) 

+ tr,^,(2y4)-4V3yf) . (Sll) 

Therefore, we obtain 

Co = 2 sgn [(y 442 _ I 2 (yf 42 j ^ (S 12 ) 

leading to the phase diagram shown in Eig. S2 (b). Eor domi¬ 
nant nearest-neighbor hybridization (phase I) we therefore ex¬ 
pect a spin texture with Wx = 1 , while we expect Wx = -1 for 
dominant next-to-nearest neighbor hybridization (phase II), 
see also Eig. 1. At the phase transition y)*^ = 2 V3y)^\ the 
hybridization vanishes along the EX line, for y)*^ = -2 V3y® 
at both the XM and XR lines, as for the r 8 model in the Main 
Text. Our numerical calculations for the r 7 model confirm 
both the MCNs and the expected spin texture, see Figs. S2(c) 
and S2(d). 


CALCULATIONS FOR THE FULL MODEL 

Assuming, as in the Main Text, that 4*^’4^^ ^ ® 

4^4 4*' ^ MCNs and therefore the spin texture depend 

on the energy of r 7 and Tg' orbitals at the X point, 
e^, respectively. If the energy difference Ae^ = ef - ef 
is much larger than the hopping between r 7 and Eg orbitals, 
Ae^ » f 7 g s max(| 481 , l 4 g 4)5 the band inversion essentially 
occurs between the (x^ - y^) and the Eg^^ orbitals. Then, hop¬ 
ping and hybridization of the E 7 orbitals is irrelevant for the 



S4 



h 


h 


(a) yW 


(b) yW 

O 




FIG. S3. Phase diagrams for the full model in the case where i/'J 
defined in Eq. (SI3) is the highest / orbital and F® = F® = 0 (a) 
and F® = F® = 0 (b). The phases I and II correspond to Wx = 1 
and wx = -1, respectively. 


depends on all four hybridization parameters, 


FIG. S2. Bandstructure without hybridization (a) and phase diagram 
(b) for the Fy model defined in Eq. (SIO) and (S8) with f® = 1, 
tf = -0.2, f® = -0.03, r® = 0.02, and ej = -3. The two dif¬ 
ferent spin textures (c-d) in phases I and II, respectively, are shown 
for hybridization parameters (F®, F®) = (0.3,0) and (F®, F®) = 
(0.1,0.1), respectively. Due to the negative eigenvalue of the re¬ 
stricted spin operator for / electrons (Eq. (14)), the spin direction 
is reversed around all HSPs when compared to Fig. 1. Note that the 
magnitude of the spin expectation value around the f is relatively 
small in (b) due to the mixing of / and d orbitals. 


Wx = sgn + 2 Vsyf ± |y'') + 3y®) 

■ (yj^’ - 2 Vsy® + 3y®)]. (Si5) 

This is a complicated phase diagram where the two different 
phases are separated by two hyper planes. For illustration pur¬ 
poses, we show the phase diagrams for the positive sign in 
Eq. (SI3), and only nearest and only next-to-nearest neighbor 
hybridizations in Figs. S3 (a) and S3 (b), respectively. 


topology and all results for the Fg model in the Main Text are 
applicable. Instead, for Ae^ «c I’ '^^e Fg ^ orbitals are 
not involved in the band inversion such that we can apply the 
results of the Fy model. 

However, at the X point there exists a non-vanishing mix¬ 
ing between those two orbitals, which we need to take into 
account for comparable energies of Fy and Fg^* orbitals. This 
leads to a continuous crossover between the two situations. 
For Ae^ = 0, we obtain an avoided crossing at the X point 
with the splitting being defined by t*g^ and Then, the high¬ 
est / band is an equal superposition of Fy and Fg*^ orbitals, 

where the sign + depends on the parameters and f^g\ 

For this orbital, the hybridization matrix with the (x^ - y^) 
at the X point to first order is given by a superposition of the 
matrices defined in Eqs. (11) and (S11); 

-i V2 d), = {(r,q,+cTyqy) + 2 V3yf ± f + 3yf) 

+ cr,q, (2y''^ - 4 VSF® + bF®) . (S14) 
Then, the winding number of the spin texture at the X points 


SPIN OPERATOR FOR / ORBITALS 


While the spin operator for the d orbitals is a block-diagonal 
matrix = 1 ® cr, this is not the case for the / orbitals [S7]. 
The spin operator for the / orbitals can be obtained by cal¬ 
culating its matrix elements for Fy and Fg states. Writing the 
spinor of / electrons as i/^ = (^/pcu /p( 2 ) , we obtain the 

spin operators 


-5 

-2V3 2V5' 

-2V3 

-9 2 vn 

2V5 

2 vn 5 


-5 

2V3 2V5 ' 

2 V3 

-9 -2 vn 

2V5 - 

-2V15 5 ^ 

-11 

1 

0 


0 

-3 0 


-4 Vs 

0 5 



(S16a) 


(S16b) 


(S16c) 


Again, the spin matrix for the reduced models is obtained by 
removing the matrix elements that include either the Fy or Fg 
orbitals. 
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